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Abstract 

Experimental realization and quantitative investigation of common-noise-induced synchroniza- 
tion of limit-cycle oscillations subject to random telegraph signals are performed using an elec- 
tronic oscillator circuit. Based on our previous formulation [K. Nagai, H. Nakao, and Y. Tsubo, 
Phys. Rev. E 71, 036217 (2005)], dynamics of the circuit is described as random phase mappings 
between two limit cycles. Lyapunov exponents characterizing the degree of synchronization are 
estimated from experimentally determined phase maps and compared with linear damping rates of 
phase differences measured directly. Noisy on-off intermittency of the phase difference as predicted 
by the theory is also confirmed experimentally. 

PACS numbers: 82.40.Bj,89.75.Da,43.50.+y 
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I. INTRODUCTION 

Synchronization of nonlinear dynanncal elements is observed in many natural systems fl 
|2j. For instance, in our body, heart cells synchronize with each other to generate heartbeats, 
and suprachiasmatic neurons synchronize with the 24 hour daily cycle to generate circadian 
rhythms 0, 4]. Many experimental investigations of synchronization have been carried out, 
e.g. using coupled chemical reactors 5), [g, Q]. Synchronization typically occurs due to 
mutual coupling or through entrainment to common periodic signals. Generally, external 
noises independently applied to the elements have negative effects on synchronization; the 
elements cannot synchronize under independent noise sources that are too extreme. 

In contrast, common or correlated external noises can synchronize uncoupled dynamical 
elements. Using neurons of rat neocortical slices, Mainen and Sejnowski [8] have shown that 
reliability of spike generation improves when a neuron receives a fluctuating input current 
compared with the case of a constant input current. This phenomenon can be considered as 
synchronization of uncoupled identical dynamical elements induced by common fluctuating 
inputs. The synchronizing effect of common fluctuating forcing, known in ecology as the 
Moran effect, describes the synchronized population dynamics of organisms due to correlated 
environmental fluctuations M. 



More explicitly, Pikovskii [lOj and Jensen [TJJ have theoretically investigated synchro- 
nization of limit-cycle oscillators induced by non-periodic external signals. Synchronization 



of uncoupled chaotic osci 
by Maritan and Banavar 



lators due to common noisy driving has been numerically studied 



12J and experimentally realized by Sanchez et al. 13j using an elec- 



tronic circuit. Synchronization (or consistency) of chaotic lasers due to common fluctuating 
signals has also been reported [14]. 

For limit-cycle oscillators, general quantitative formulations of common-noise-induced 
synchronization can be developed using the phase-reduction method [lj. Teramae and 
Tanaka [l5j have proven that the synchronized state of uncoupled limit cycle oscillators 
subject to a common weak Gaussian noise is always statistically stabilized, and their theory 
has been further gen eralized to provide global stability of phase coherent states induced by 
correlated noises 16]. The cases where limit-cycle oscillators are stimulated by a common 

rr n 

telegraph noise [17] or by a common impulsive noise [18( have also been investigated theo- 
retically using random phase-map descriptions. Recently, synchronization due to common 
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random impulses has been studied experimentally with an electronic circuit, and some of 
the theoretical predictions have been quantitatively verified jl^ . 

In this paper, we experimentally investigate common-noise-induced synchronization using 
an electronic circuit undergoing periodic oscillations. As the random signal, we use a random 
telegraph noise, which is the simplest example of colored non-Gaussian noises; it can easily 
be generated in experiments and facilitates analytical treatments. In this case, we have two 
limit cycles corresponding to two values of the driving signal, in contrast to the previous 
experiment using random impulses where the system possessed only a single limit-cycle 
orbit 19]. When the switching interval of the driving signal is sufficiently long, the circuit 
state is mostly on either of the limit cycles at the moments of switching, so that its dynamics 
can be described in terms of phase mappings between the two limit cycles. We experimentally 
determine the phase maps of the electronic circuit and estimate the Lyapunov exponents 
characterizing the degree of synchronization from the phase maps based on our previous 
theory. The Lyapunov exponents are then quantitatively compared with the damping rates 
of small phase differences measured directly. We also confirm that noisy on-off intermittency 
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re phase difference, which is typically expected for random- mapping systems 13, 



2d, 



21 



231 ]. actually occurs in our electronic circuit. 



II. EXPERIMENTS 



A. Setup 

The experiments were performed using an electronic circuit shown schematically in 
Fig. [T](a), where an LM741 was used as the Op- Amp and the circuit parameters were set as 
follows: Ri = 100 kQ±5 kfi, R 2 = 81 kfi±4.5 kfi, R 3 = 470 kfi±23.5 kfi, R 4 = 100 Q±5 Q, 
C\ = 100 nF ± 5 nF, and Ci = 10 /xF ± 1 /xF. Voltages of positive and negative power sup- 
plies to the Op-Amp were fixed at ±3.0 V (0 V indicates the ground voltage) using a DC 
power source (PMM18-2.5DU, Kikusui Electronics Co.). Note that we use the Op-Amp 
under positive feedback conditions to generate oscillations, so that the golden rule of an 
Op- Amp (V+ = V-) does not hold in our experiments. 

The source voltage of the MOSFET (2SK2201, Toshiba Co.) was fixed at -4.0 V with 
another DC power source (E3630A, Hewlett Packard) and the gate voltage V g (t) was con- 



3 



trolled by the external signal. Voltage traces V+(t) and VI (t) were measured from the 
circuit as shown in Fig. Ufa). Control of V g (t) and measurements of V+(t) and V-(t) were 
performed with an AD/DA converter AIO-163202F-PE (Contec Co.). When V g (t) was fixed 
at a constant value between -6.0 V and -2.45 V, the circuit exhibited limit-cycle oscillations. 
Figure [T](b) shows a limit-cycle orbit on (V+, VJ) plane at V^(£) = —6.0 V, and Figs. [CJb) 
and (c) display the corresponding time series of V+(t) and VL(t). 



(a) (b) 




FIG. 1: (Color online) Experimental setup, (a) An electronic circuit used in the experiment, (b) A 
limit-cycle orbit observed at V g {t) = —6.0 V on the (V+, V_) plane. Average period was 0.0445 s. 
(c) Time series of V+(t) (left) and V-(t) (right). 

We repeatedly switched the gate voltage V g (t) between two values V g \ and V g2 to simulate 
a random telegraph signal. The switching events obeyed a computer-generated Poisson 
process, namely, the switching interval D was an exponentially-distributed random variable 
with mean interval r = 0.2 s. We generated D by the formula D = — l/rlog(l — u) 
with u being computer-generated pseudo random numbers in [0,1]. The typical relaxation 
time of the circuit to converge to either of the limit cycles was shorter than 0.01 s. By 
applying the same time sequence of V g (t) to the circuit repeatedly, we performed consecutive 
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measurements of the time series of V+(t), and then repeated this procedure with different 
realizations of the random telegraph signals. In the experiment, V g i was fixed to -6.0 V and 
V g 2 was varied between -6.0 V and -2.45 V. 

We defined a phase <p(t) of the circuit state from the time series of V + (t) as follows. The 
origin of the phase (0 = 0) was taken as the moment when V + (t) changed its sign from 
negative to positive. Each time V+(t) crossed V from negative to positive, the phase was 
reset to 0. Between successive zero-crossing events, the phase was increased with a constant 
frequency from to 1. Note that the frequency of the oscillator was not constant but changed 
from cycle to cycle due to random switching of the driving signal. From the two consecutive 
time series of V + (t), we obtained time series of the absolute phase differences \A<f)(t) | between 
two experimental trials (restricted to [0, 1] using the periodicity of the oscillator) under the 
same time sequence of V g (t). 

B. Results 

Figure [2](a) shows the time evolution of the absolute phase difference |A0(t)| between two 
experimental trials under a constant input, V g (t) = V g \ = —6.0 V, which increased linearly 
with time t. Even after fine tuning, we observe that the average period of oscillations differed 
slightly across experimental trials (±0.1 %). 

Figures E^b-l) and (b-2) show the time series of A0(t) observed under the common 
telegraph noise, where V g (t) was switched between two values (V^i = —6.0 V and V g2 = 
—2.5 V). Large changes of |A0(i)| were observed only in short time windows after the 
switching events of V g (t) as shown in Fig. [2](b-l). The two experimental trials driven with 
the same input signal became mostly synchronized, but occasionally there were interruptions 
by short desynchronization events as shown in Fig. [2](b-2). 

To characterize this characteristic intermittent behavior of the phase difference, we mea- 
sured the distribution of laminar intervals during which |A0(i)| was smaller than a certain 
threshold, |A0(t)| < 0.01. As shown in Fig. [2]^c-l), the stationary distribution of laminar 
intervals appeared to follow a power law, whose exponent was approximately —1.5. We also 
measured the stationary distribution of the absolute phase difference |A0(t)|, which also 
exhibited a power-law tail with an exponent roughly —1 (Fig. [2]^c-2)). Though the phase 
difference |A0(t)| shown in the figure is restricted to the range [0,1], the phase difference 
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occasionally exhibited jumps of magnitude 1 due to phase slippage of one period (e.g. near 
t = 1200 and t = 1600 in Fig. |2j(b-2)) [2411 . 



Thus, when the gate voltage V g (t) was switched between two values randomly, different 
experimental trials tended to be synchronized even under the effect of slight differences in 
average periods and experimental noise. As we explain later, the characteristic behavior of 
the phase difference |A0(t)| was due to noisy on-off intermittency 
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< 6i(V+, VJ) < 1 for LCI and also in the 
frequency 1/Ti with Ti the period of LCI 



0,1s 



III. ANALYSIS 
A. Theory 

Here we briefly summarize our previous theory on the synchronization of uncoupled os- 
cillators driven by a common random telegraph signal [l?]]. Corresponding to the two dif- 
ferent values of the input signal V g (t), the circuit exhibits two different limit cycles, LCI 
for V g {t) = V g i and LC2 for V g {t) = V g 2. When V g (t) = V g ±, we can define a phase 

[V, VJ) plane, which increases with a constant 
3j. Similarly, when V g (t) = V g 2, another phase 
< 9 2 (V + ,VJ) < 1 can be defined in the (V+,V_) plane that increases with a constant 
frequency 1/T 2 , where T 2 is the period of LC2. The origins of 9i and 82 are taken as the 
points where V + (t) crosses V from negative to positive on LCI or LC2. Combining these, 
we introduce a new phase 9{V + , V_, V g ) of the circuit state as 

x f 9i(V+, VJ) (when V g = Ki), 
9(V + ,V„,V g ) = I V + ' K g g ; ' (1) 

( 2 (V+,VJ) (when V g = V g2 ). 

Note that this 6 is different from the phase <fi that we defined in the previous section 
by linearly interpolating successive zero-crossing events. 9 jumps discontinously at the mo- 
ments when V g (t) switches, because 9\ and 82 increase with strictly constant frequencies. In 
contrast, is continuous even when V g (t) is fluctuating (except the zero-crossing events of 
V+(t)), but its frequency differs from cycle to cycle. When V g (t) is kept constant for longer 
than one period of oscillation and than the relaxation time of the circuit to the limit cycle, 
coincides with 9. This difference in the definition of the phase variables yields only small 
bounded discrepancies in measuring the phase differences between two time series. 
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FIG. 2: (Color online) (a) Evolution of the absolute phase difference |A</>(i)| under a constant 
external signal V g (t) = —6.0 V depicted using doubly logarithmic scales. The solid curve represents 
the experimental data and the broken line with unit slope represents linear dependence on t. (b) 
Time series of the phase difference, |A</>(i)|, with V g (t) a random telegraph signal (V s i = —6.0 V 
and V g 2 = —2.5 V). (b-1) Short time series of V g (t) (top) and |A0(i)| (bottom), (b-2) Long time 
series of |A^>(t)|, exhibiting noisy on-off intermittency. (c-1) Distribution of the laminar length. 
The broken line represents i -1,5 . (c-2) Distribution of the absolute phase difference |Ac/>(i)|. The 
broken line represents |A^>(t)| _1 . 
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FIG. 3: (Color online) Phase mappings between two limit cycles, (a) When V g (t) switches from 
V g i = —6.0 V to V g 2 = —2.45 V, points on LCI (solid curve) for Vgi are mapped to the corre- 
sponding points on LC2 (broken curve) for V g 2 as indicated by arrows, (b) From V g 2 = —2.45 V 
to V gl = -6.0 V. 

We assume that the average switching time of the input signal, r, is sufficiently longer 
than the relaxation time of the orbit to LCI or LC2 at fixed V g (t). In our experiments, 
we took r = 0.2 s and the relaxation time was typically shorter than 0.01 s, so that this 
condition was satisfied. The orbit of the circuit is then almost always on one of the limit 
cycles when V g (t) is switched between two values. Therefore, we can describe the dynamics 
of the circuit under randomly switched V g (t) as alternating phase mappings between 6\ and 
9 2 as shown in Fig. [3j We denote the mapping from #i to 62 that takes place when V g (t) 
switches from V^i to V g 2 as 6 2 = f 12(61), and the mapping from 62 to Q\ when V g {t) switches 
from V g2 to V gl as 9 X = / 21 (0 2 ) [17(] . 

Let us denote the phase on LCI just before V g {t) switches from V g i to V g2 for the n-th 
time as and the phase on LC2 just before V g (t) switches from V g 2 to V^i for the n-th time 
as 6V;. Then the phase dynamics of the orbit can be described as (assuming V g (t = 0) = V^i) 

^ +1 = /2i(^ n ) + ^ +1 , ^ = /i 2 W) + 4, (2) 
where s™ and s 2 are exponentially distributed random switching intervals whose probability 
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distributions are given by 

= ^exp (-^) , P 2 (4) = ^exp , (3) 

respectively. Small deviations A8™, A9^ from 0™, 9% obey linearized equations 

A^ +1 = f'M)MZ> A6% = f[MWi, (4) 

where /' denotes the derivative function of /. After the n-th switching, the amplitude of 
the deviation A6™ is given by 

|A*?| = (j[ \fM)\ J |Ae}|, (5) 

and |A^2 I similarly. Thus, for large n, 

|A0?| |A02? 



|A6>J| - |A^" 6Xp(An) (6) 
holds, where the Lyapunov exponent A is given by 

A= lim -lnfn \f'M)\ \fM)\ ) =Ai + A 2 (7) 

with 

Ai~ /V In |/{ 2 (0i)|, A 2 ~ fd9 2 In i/^(0 2 )|. (8) 
Jo Jo 

Here, we approximate the average over the stationary distribution of the phase 9\ or 62 under 

the effect of telegraph noises by the average over the uniform distribution in each equation, 

because the phases 9\ and 6* 2 are almost uniformly distributed on LCI and LC2 when r is 

sufficiently larger than T\ and T 2 [171 ] . In our experiments, we used r = 0.2, whereas the 

period of oscillations was about 0.05 s. Therefore, this condition was satisfied. 

When A = Ai + A 2 < 0, phase synchronization induced by the random telegraph signal 

is expected. Note that the switching step n is approximately related to the real time t as 

n ~ t/2r, so that 

|A9(i)l -expfAA (9) 



|A9(0)| r \2t 

holds for large n, t. 
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B. Determination of phase maps 



We experimentally determined the phase maps 9\ = f 21(02) and 9 2 = /i2(#i) as follows. 
We first measured the period 7\ of LCI under a constant input signal V g (t) = V g \. After 
setting V g (t) at V g i and relaxing the circuit for 0.5 s, intervals between successive zero- 
crossing events of V + (t) from negative to positive values were measured for 5 s. T\ was 
determined by averaging these intervals. We then measured the period T 2 of LC2 at V g (t) = 
V g 2 in a similar way. 

When the measurement of T 2 was completed (this moment was defined as t = 0), the 
following regular telegraph signal (shown schematically in Fig. H]) was applied as the probing 
input: 



Vgi (t*- 1 < t < t\) , 



V g (t) = { 



(10) 



v s2 {t\ < t < 4) , 

for i = 1, • • • , 200, where the i-th switching time from Vgi to V g 2 is given by 

t\ = d 1 , 4 = ^ + ^24 (i>2), (11) 

k=l 

and the subsequent i-th switching time from V g 2 to Vgi is given by 

i 

4 = ^24 (z>l). (12) 

k=l 

Here, di denotes the length of the i-th constant interval of the input signal. To avoid 
undesirable synchronization with the probe signal, we gradually increased di as di = 0.5 + 
0.001(i - 1) (1 < i < 200). 

Using this V g (t), we measured the phase < <I?\ < 1 of the circuit at t\, 

{(<2>1,4) | i = I,--- ,200}, (13) 
and the phase < <P 2 < 1 of the circuit at tf 2 , 

{(^ 2 ,4) |i = l,-..,200}. (14) 
These phases <P[ and <P 2 can be identified with the phases 6\ and 62 used in the theory via 

= 02(t = t 2 )=<P 2 . (15) 
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FIG. 4: (Color online) The probe signal used to determine the phase maps. Before t = 0, T\ and 
T 2 were measured. After t = 0, the phases were measured at t| and t\ for 1 < i < 200. 

Note that the two types of the phases coincide here because the constant intervals di of the 
probe signal are always longer than T 1; T 2 , and the relaxation time of the circuit to LCI or 
LC2. 

At each switching event, the phase jumped from <P\ to <P l 2 (V g i — * or from <f>\ to <I?\ 
(V g 2 — > V g i). The destination phases <P\ and <P 2 just after t\ and t 2 were obtained from <P\ 
and <P 2 as 



si+l 



) mod 1, 



— ^- ) mod 1, 

To 



(16) 



where the definitions were made modulo 1 to restrict the phases to [0, 1]. Thus, we obtained 
200 realizations of the phase mappings, 



{$*-<Pi|i = l,-..,200}, {$i-<^ |i = l,-.. ,200}. (17) 
We constructed the raw phase maps by piecewise-linearly interpolating these data as 



(18) 



which were still non-smooth functions due to experimental fluctuations. 

Generally, the phase map has a trivial diagonal component, namely, the identity- map 
component that exists even when V s i = V g 2, and additional non-trivial components reflecting 
the nonlinear transition dynamics between the limit cycles. We estimated the underlying 
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smooth phase maps fu and /21 from the raw phase maps /[| w and f^f by low-pass filtering 
using the 20 lowest Fourier modes as 

20 

/ia(0i)=0i+ 4exp(2 7 r^ 1 ), 

fc=-20 
20 

/ 21 (0 2 ) = # 2 + ^ 4 exp(27rA;^ 2 ), (19) 

fc=-20 

where cj 2 and 4 are Fourier coefficients of the non-trivial components of /{| w and /|i w , 
defined as 

00 

AT(^i)-^i= Yl 4exp(2 7 r/c^ 1 ), 
/c=— 00 
00 

/ 2 T(^ 2 ) - <2>2 = ]T 4 exp(2vr^ 2 ). (20) 

fe=— 00 

Figure [5](a) displays examples of the phase maps obtained with the above procedure. 
C. Lyapunov exponents and damping rates 

From the experimentally determined phase maps /i 2 (#i) and / 2 i(6' 2 ), the Lyapunov ex- 
ponent A = Ai + A 2 can be estimated. This A can be compared with the damping rate rd a m P 
of small phase differences between two trials subject to the same telegraph noise. As we 
have already explained, the difference in the definition of two phases results in only a small 
bounded discrepancy between A8(t) and A(f>{t), so that it does not affect the Lyapunov 
exponents or the damping rates. 

When the phase difference A0(t) is small, we expect the ensemble average of A0(t) over 
many realizations of the random telegraph signal to shrink exponentially as 

(|A0(t)|) ~exp(r damp t), (21) 

where rdamp < is the damping rate. As the average switching interval of V g (t) is r, 

'"damp = 7; (22) 

zr 

will approximately hold for large t and n provided that the previous analysis based on the 
phase-mapping description is reasonable. 
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FIG. 5: (Color online) (a) Experimentally determined phase maps. Circles show the raw data 
($\, <P\ — $\) (left) or (# 2 , $\ — ^2) ( r ight), and curves represent the non-trivial part of the estimated 
phase maps fui^i) — Oi (left) or /2i(#2) — &2 (right) obtained after low-pass filtering. V g \ = —6.0 V 
and V g 2 = —2.5 V. (b) Comparison of the Lyapunov exponents A = A1+A2 with the linear damping 
rates r^mp directly measured from \Acj>(t)\. The line represents 0.4 r^mp = Ai + A2. Filled circles 
represent the data obtained by varying V g 2 (—2.55 V < V g 2 < —2.45 V with intervals of 0.01 V) 
while fixing V„\ at —6.0 V. 
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Figure [5(b) compares the Lyapunov exponent A with the damping rate rdamp obtained 
for different values of V g 2 with fixed V E \. The damping rate rdamp was directly measured 
from |A0(t)| as 

1 

rdamP " M-E^AT fc 
where A0 upper = 0.2, A0 lower = 0.05, AT k is the time needed for |A0| to be damped from 
A0 upper to A0i ower , and M is the number of such shrinkage events in the time series of \ A<f>(t) \ . 
We measured 20 time sequences of V+(t) for 120 s and calculated 19 time sequences of the 
phase difference A<fi(t) between two consecutive time sequences of V+(t) to obtain r damp . As 
shown in Fig. [5](c), pairs of (A, rdamp) estimated for various values of V g 2 approximately fall 
on the straight line A = 0.4 rd am p ( r = 0.2 in the experiment), which quantitatively verifies 
the validity of the phase-mapping description of our experiments. 



A0; 



lower 



upper 



(23) 



D. Noisy on-off intermittency 



We have focused so far on the average behavior of the phase difference. The phase dif- 
ference |A0(t)| decreases on average when A = Ai + A2 < 0. However, as can be seen from 
Eq. (jSJ), small phase differences A9i and A62 are driven multiplicatively by the random ap- 
plication of two phase maps. This is a typical situation where noisy on-off (or modulational) 



intermittency is expected over long time scales [17|, [20j, [21|, |22|, [23j; due to small noises or 
heterogeneity inherent in the system, individual time sequence of |A0(t)| can occasionally 
grow due to random multiplication even if A = Ai + A2 < 0, resulting in repetitive transient 
bursting. As already shown in Fig. 12(c), this is the case for our electronic circuit. The power- 
law distribution of the laminar interval with the exponent —1.5 as shown in Fig. |2(c-l), and 
the power-law distribution of the amplitude of the phase difference as shown in Fi g. BT c-2) 



are consistent with the theoretical predictions on noisy on-off intermittency 



21 



22 



23|. 



IV. CONCLUSIONS 



We investigated synchronization between different experimental trials induced by common 
telegraph noises using an electronic circuit undergoing limit-cycle oscillations. The dynamics 
of the circuit could be described in terms of random phase mappings. We experimentally 
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determined the phase maps and quantitatively verified that the Lyapunov exponents deter- 
mined from the phase maps agreed with the damping rates measured directly from the time 
series of small phase differences. We also confirmed that noisy on-off intermittency of the 
phase difference actually occurs. 

The mechanism leading to synchronization that we demonstrated using an electronic 
circuit in this paper is general and is expected to be observed in various systems undergoing 
limit-cycle oscillations. 
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